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Abstract 

Anderson transition in three-dimensional systems is investigated using 
renormalization group theory. /3-function of a very simple form is derived 
from a self-consistent consideration, and it gives a value 1 + 1/V3 = 1-58 
for the critical index of the transition, which is very close to those obtained 
by numerical studies. 
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1 Introduction 

Metal-insulator transition of electronic systems due to randomness was pre- 
dicted by Anderson in 1958. jlj Let i be a quantity which controls the conduc- 
tivity, such as impurity density or electron density. Then the system undergoes 
a transition from a conductor to an insulator at x — x c , the critical value of x 
(in the following we assume that the system is metallic for x > x c ). One of the 
main interests on this subject is the critical behavior of the transition, namely, 
how the conductivity cr{x) depends on x near the critical point. In spite of an 
amount of theoretical studies, little progress has been made about this problem, 
until the appearance of the renormalization group (or scaling) theory developed 
by Wegner M and by Abrahams et al. f| . 

The scaling hypothesis by Abrahams et al. leads to the renormalization 
group equation |J] 
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where g(L) is the conductance of a system of size L. The renormalization 
group theory predicts that the transition is continuous in three-dimension: u(x) 
vanishes continuously at x — x c . The critical index (or exponent) v of the 
transition is defined by the equation 

a(x) = C(x - x c ) v . (2) 

Here C is a constant which is dependent on the microscopic structure of the 
system, while v is believed to be universal and dependent only on the univer- 
sality class, i.e., the fundamental symmetry of the hamiltonian. Below we will 
treat the orthogonal class, in which the Hamiltonian is composed of the kinetic 
energy and the impurity potential, without spin-orbit coupling nor magnetic 
field. 

The critical index v is given in terms of the (3- function (3(g): 

where the prime on (3(g) indicates the derivative, and g c is the fixed point value 
of g, satisfying 

(3(g c ) = o. (4) 

Elaborate calculations of the /3-function have been done by Wegner || and 
Hikami, || and their results lead to v < 1. This problem was investigated 
also by numerical methods. Kramer and macKinnon developed a method to 
estimate v from numerical data, and they obtained a value considerably larger 
than unity, J?], |) Calculations with high accuracy have been performed recently 
by Slevin and Ohtsuki, and they found that J|, [l(J 

i/ = 1.57 ~ 1.58. (5) 

Thus the discrepancy between the analytic and the numerical methods is 
significant. At present, as long as the accuracy of the results are concerned, 
the analytic methods seem less reliable than the numerical methods. However, 
it is still important to pursue better analytic methods, because there are some 
aspects of which the better understandings are obtained only by analytic meth- 
ods. 

In this paper we propose an analytic theory based on an idea different from 
those of existing theories. We will find it well reproduces the results of numerical 
studies. 

2 Perturbational Calculation of /3-function 

In this section we briefly review the method applied so far to calculate (3- 
function, for it is the starting point of the present theory, too. The details 
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of the calculations are reviewed in Ref. j| To be specific, we consider an elec- 
tron gas system with randomly distributed impurities. The hamiltonian is of 
the form 
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^i^+E^-^)- ( 6 ) 

i 

where RiS are the positions of the impurities, which arc assume to be dis- 
tributed uniformly with the average density Cj, independently of each other. 

So far the /3-function was calculated in the form of the power series in l/g. 
The term of 1 /<? arises from the lowest order interference effect of electron wave 
scattered by impurities. Without the interference correction, the conductivity 
is given by Boltzmann (classical) conductivity 

_ n e e 2 T 

a = <j B = , (7) 

m 

where n e is the electron density, and r is the electron scattering time by impu- 
rities given by 

- = -r-c iU 2 N 8 
r n 

No being the density of states at the Fermi level. 

In thermal Green's function formalism, [[llj the interference correction to the 
Boltzmann conductivity comes from the term corresponding to the Feymann 
graph shown in Fig. |l[ Here the solid lines and the dotted lines indicate the one 
electron Green's functions and the impurity potential, respectively; the wavy 
line indicates the cooperon propagator. Up to this order, the conductivity of an 
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Figure 1: The Feynmann graph for the lowest order interference correction 
infinitely large system is given by 

dq 
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where Do is the diffusion coefficient defined by 
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vf being the Fermi velocity, and the integral is to be done in the region q < 
q c S3 1/vft. 

The conductance g(L) is defined in terms of the conductivity cr(L) of the 
system of size L as 



j(L) = rLa(L) , 



(11) 



where r is a constant of dimension of resistance to make g(L) dimensionless, to 
be determined later. We postulate that <x(L) is given by the right hand side of 
eq. (0) in which the integral is cut off at q — a/L, a being a numerical constant 
of order 1 (we consider the region a/L < q c ). Then we obtain 



g(L) = rL 



e 2 / a\ 
aB -—h\ qc -L) 



and /3-function is obtained from eq. (hi): 
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(3(g) = i 



(13) 



where we have taken r so that 
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As is seen from eqs. (||) and (||), the ambiguity of r within a numerical factor 
does not affect the value of v, and we find that 
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Calculations of /3-function up to the order 1/g 4 were done by Wegner || 
and Hikami, ||. Their /3-functions lead to v = 0.56, and the deviation from 
the numerical value becomes even larger. The inclusion of the higher order 
seems very difficult. Moreover, it is not certain if the 1/g expansion is a proper 
expansion. Thus we have to investigate the problem from a different point of 



3 Self-Consistent Treatment 

Self-consistent treatments of Anderson transition were developed by Vollhardt 
and Woelfle, (T^| and by the present author. jl3| Although those theories con- 
tributed to the deeper understanding of the problem, they give v = 1 like 



the renormalization group theory with the /3-function (13). Those theories are 



based on the requirement of the consistency between a and Dq in the integrand 
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(cooperon propagator) in cq. (||) for L — ► oo. In this paper we require the 
consistency including the L dependence of these quantities. 

In this sense eq. (|]) is not self-consistent, for a{L) depends on L while Dq is 
a constant. The diffusion coefficient must be independent of q for small q, but it 
can depend on q for large q. Therefore, we assume that the diffusion coefficient 
behave like oc q 1 for large q. We consider the case when a/L is large in some 
sense. Then, assuming that the cooperon propagator is proportional to l/g 2 " 1 " 7 , 
and putting it into eq. (|^) we find that 



(16) 



if we neglect the terms independent of L. We identify a/L with q, and from 
the self-consistency requirement we find that 7 = 1/2 and that the diffusion 
coefficient should behave like ^fq for large q. Scaling arguments predict a dif- 
ferent q dependence for large q, fLlfl and this problem will be discussed in the 
last section. 

Thus, among various possibilities, we assume a simple form for the cooperon 
propagator: 
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D q 2 (l 

where c is a numerical constant to be determined self-consistently, and 

A„ 1 
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The combination of q with Ao is based on the scaling assumption that a length 
should be scaled by A . 

Then, by replacing the integrand in eq. @ with the expression (fl7|), we 
obtain 



D(L) 
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log 
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Cy/a\ Q /L + 1 



(19) 



In comparing it with the denominator of the expression (^7|), we neglect 
the last two terms in the square bracket. In fact, we will see that the terms 
independent of L does not affect the renormalization group equation. Moreover, 
we have to compare them in the regions a\o/ L > 1 and the L dependence 
of the logarithmic term is weaker than that of \J aAo / L. Thus, identifying a/L 
with q, the self-consistency requirement gives 
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4 Renormalization Group Equation 



We will derive a renormalization group equation from eq. (|l9|). Here we define 

L 



9{L) = 
X(L) = 



A(L)' 
1 



2whN D(L) ' 



(21) 
(22) 



where 6 is a numerical constant, and we easily find that this is equivalent to 
eq. (Thy). Then we obtain 



L 



V2 



g(L) =6—1 VXoq c 



V2 



b \ a- h 6— log — ; — , 

TT V A A v/aAo/L + tt/V2 



and 



1 



dlogg(£) _ 

dlogi gr(i) 



9{L) ~V\ b— 



Ao Ao 2 + rj- 1 y/bL/X 



(23) 



(24) 



with r\ = \/ba I {\[2ti) . 

We assume that the Anderson transition is described by renormalization 
group equation. Then, the right hand side of this equation have to be a universal 
function only of g(L). From eqs. (|l8|), (|2|) and (|22|), we find that bL/Xo is <?(£) 
without the interference correction. Therefore, it is reasonable to identify bL/Xo 
with exact g(L). With this replacement, the /3-function is obtained: 



^ = 1-^ + ^-7= 



(25) 



where we have chosen 6 so that rj = 1. 



5 Critical Index 

The critical index ^ is obtained from eqs. (||) and (Q). Using eq. ( p5| ) we easily 
find that 



, 9c = 4-2V3, 

f = 1 + -= = 1.58. 

V3 



(26) 
(27) 



This value of ^ agrees with that of numerical studies, || [T(| i.e., v = 1.57 ~ 1.58. 
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6 Summaries and Discussions 



On the basis of a self-consistent consideration, we have derived a renormalization 
group equation for Anderson transition. It gives a critical index very close to 
those by numerical studies. We have not intended to derive an exact result, and 
such a good agreement is not the most important result of the present theory. 
It should be noted, however, that so far the analytic methods have never given a 
critical index comparable with those by numerical methods, within the author's 
knowledge, except for the rather phenomenological theory by Shapiro. jl5L |Tg| 

As regards the one parameter scaling, the recent numerical study by Ohtsuki 
and Slevin revealed that the critical index is little dependent on the definition of 
g(L), i.e., what kind of averaged conductance we take for it. JItJ This result is 
very important because it indicates that a simple one parameter renormalization 
group theory is reliable. 

A crucial assumption in this theory is the form of the cooperon propagator 
given by eq. (|l7|). As was mentioned in Sec. [| scaling arguments predict 
that the diffusion coefficient D(q) should depend on q like D(q) oc q 1+r) with 
rf w 1.3 at the critical point, or for such q that qX 3> 1, where A is the coherence 
length defined by eq. (p2|). On the other hand, from eq. ( |2l| ) we find that 
g w L/X « l/(gA), and as regards the critical behavior of the transition, the 
relevant values of g are g c ~ 1- Therefore, the relevant region of q is qX w 1, 
and the assumption D(q) oc yfq in this region the does not contradict to the 
scaling arguments. 

Since the /3-function ( p5| ) can not be obtained by the expansion in l/g, 
probably it is impossible to derive the ^fq term in the expression (^) by per- 
turbational approaches. The f3- function, however, reduces to eq. ([O) for g ^> 1, 
and it might be the reason why the 1/g expansion seems to work for first order. 
For small g, on the other hand, the /3-function ( |25| ) does not reduce to an ex- 
pected form (3{g) log g, but, it does not invalidate the present results, for, as 
was mentioned above, the important region of g is g « 1. 

In this context, the weak localization theory, which was very successful in 
explaining the magnetoresistance in fully metallic regions j|, |l8[ |2(|, is not 
affected by the introduction of ^fq term into the cooperon propagator, for the 
relevant scale of q in that theory is 1/£b — y eB/h, B being the magnetic flux 
density, and is much smaller than 1/Aq. 

As for the experiments, many of them suggest that v w 1. ]2l| However, 



recent careful analyses by Itoh et al. indicate a possibility of v > 1. []22| From a 
theoretical point of view, it is important to clarify the roles of electron-electron 
interaction. 
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